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ABSTRACT 


Instability in liquid monopropellant rocket motors may be 
corrected by the use of a feedback servomechanism. This 
mechanism consists essentially of a pressure pickup which senses 
pressure oscillations in the combustion chamber, an amplifier and 
a variable capacitance in the feed line. It is shown that a feedback 
system with an apbittary sampling circuit which causes the capacie 
tance in the line to complete its ‘ot cycle of variation once for 
every several cycles of combustion pressure oscillations can be 
made to stabilize the oscillations for all values of combustion time 
lay for a particular motor. It is believed that this system of 


stabilization may be applied to monopropellant motors in general. 
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SYMBOLS AND DEFINITIONS 


t = time 
cc = time lag 
p = instantaneous value cf pressure 


= dimensionless quantity relating pressure and time lag 


P = steady-state pressure in the combustion chamber 

Ss = injector pressure drop 

EE = p/ZOp , injector pressure drop parameter 

f = distance along feed line between pump outlet and injector 
m = stendy-state mass flow rate 


= cross-sectional area. of the feed line 


gas residence time in the combustion chamber 


“4H OD » 
Lo 
| 


= 0m /ZO6 Aq line inertia characteristic parameter 


instantaneous value of pump outlet pressure 


So 
\e) 
it 


steady-state value of purnp outlet pressure 


OG} 
Oo 
" 


= instantaneous value of pump outlet mass flow rate 


- rm ( p,- 6) i e, (msm) pump delivery characteristic parameter 


& Q 5: 
it 


= (p-e )/P » fractional variation of pressure in the combustion 
chamber 


t/84 , reduced time 


‘ae 
it 


2 = t/Qg » reduced time lag 


iM “2 (m;—m )/en fractional variation in injector mass flow rate 
~ = reduced amplification coefficient 
Ww = reduced angular frequency (real frequency divided by resi- 


dence time) 
= = N+ ow 
*.(%) s instantaneous value of capacitance 


ce ) 


transfer function of sampling circuit 


tt 


pee transfer function of amplifier 


I. INTRODUCTION 


The low frequency, ''chugging", oscillations in liquid rocket 
motors have been attributed to the coupled effects of the feed system 
and combustion chamber dynamics. According to D. F. Gunder and 
D. R. Friant {}), the pressure oscillations can take place in the 
combustion chamber merely as a result of small changes in the 
velocity of flow in the injected fuel; that is, if for any reason the flow 
rate decreases, the rate of burning and hence the pressure will 
decrease a short time later in the chamber. Then, as the pressure 
drop is felt at the injector the flow rate will increase and a short time 
later the pressure at the point of combustion will increase resulting 
finally in a decrease in rate of flow from the injection nozzle, the 
cycle continuing indefinitely. The explanation for this phenomenon 
involves the concepts of time lag between the instant of injection and 
transformation of the fuel to hot combustion products, and of the 
dependence of the rate of burning on combustion chamber pressure. 

M. Summerfield (2) has shown how the length of feed line, the 
velocity of the propellant in the feed line, the ratio of feed pressure 
to chamber pressure and the ratio of chamber volume to nozzle area 
affect chugging oscillations and has discussed the trend of changes 
required in the feed systern and motor parameters to suppress these 
oscillations. 

or Cxteecle? has introduced the concept of variable time lag 
basing this concept on a qualitative analysis of the complicated 


processes that affect the propellants from the time of injection to the 





time of burning. He has formulated a relation between time lag 

and pressure showing that if the pressure is oscillating, the time 
lag is also an oscillating quantity around an average value. Using 
this idea he has demonstrated that self-excited oscillations in the 
combustion chamber can exist even in the absence of any variation 
in the injection rate. In the case where the time lag is decreasing 
the burning of the particles that were injected later catch up with the 
burning of those that were injected earlier and an increase in rate of 
burning would result, the opposite being true if the time lag were 
increasing; then, if the variations in time lag coincided with the 
variations in pressure in a certain manner, self-excited oscillations 
would be produced. 

H. S. Tsien (4) has demonstrated that a feedback servomech- 
anism may be used to stabilize combustion chamber pressure 
oscillations for all values of time lag, thus making it possible to 
stabilize motors for which a change in configuration parameters alone 
would not correct instability or where such changes would not be 
feasible. The servomechanism Tsien proposed consists essentially 
of a feedback servocontrol which senses the oscillating pressure in 
the combustion chamber and changes the propellant flow by the proper 
amount and in the proper phase to damp the oscillations. Sucha 
system is shown schematically on Figure l. 

A servomechanism in which the variable capacitance is 
required to oscillate at the same frequency as the chamber pressure 


fluctuations is likely to produce at least two important problems 
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arising from the rate at which these oscillations generally occur, 
namely in the order of 100 per second. One is the difficulty of pro- 
ducing sizeable changes in capacitance at such frequencies, and the 
other is cavitation in the fluid. If the variable capacitance consists 
of a piston-cylinder mechanism, it would help to lessen the problems 
to move the piston down on the compression stroke at the same 
frequency as the chamber oscillation and to withdraw the piston 
slowly to the maximum-capacity position in the time required for 
several pressure oscillations in the chamber to take place. With the 
proper phasing and amplitude it should be possible to damp out 
pressure oscillations in the chamber. 

The purpose of this investigation is to design a feedback 
circuit that will detect symmetric (say sinusoidal) oscillations in the 
combustion chamber and provide signals to the servomotor that will 
produce the asymmetric motion described above and provide stability 
for all values of combustion time lag. 

As an extension of Tsien's work, F. E. Marble and 
BD. W. Cox (9) have shown that a feedback servomechanism mey 
stabilize liquid bipropellant motors as well as monopropellant motors. 
Their general method for designing the stabilizing transfer function 


of the feedback loop will be used in this investigation. 





li, DYNAMICS OF THE MONOPROPELLANT MOTOR 


The basic motor without feedback control considered in this 
investication essentially consists of a combustion chamber, injector, 
feed line and fuel pump. it may be thought of as two dynamic systems: 
one, the combustion chamber and the discharge nozzle; the other, the 
propellant supply system consisting of the propellant purnp, propellant 
line and injector. These two systems are coupled through the require- 
ment that the mass of propellant discharged from the propellant line 
appear in the combustion chamber. The performance of the combustion 
chamber is characterized by a mean pressure p , a mean residence 


time 9, anda mean flow rate ro. When the chamber pressure 
P- Pp 
p 
The performance of the propellant feed system depends upon the cross- 


varies from p its fractional variation will be denoted = ; 
sections] area A andlengsth ? of the feed line, the pressure drop Aj 


across the injector and the local slope a oi the mass flow-pressure 


— 


characteristics of the pump where a= -— ae eee —Re-_ ie 
ae Po Ps 
fractional variation in propellant injection rate will be defined as 
ae 2 
M = —_! 
a) 


The fuel injected into the chamber does not burn immediately 
but only after it has remained a certain time t inthe chamber. The 
time «+ required for the transformation of propellant into products 
of combustion rnay be represented by means of the integral 

it 
+ (6, a = constant 


t-t 
where f(p, T A is a function depending upon the mechanics of the 


process, and Tyg is the ambient gas temperature in the neighborhood 
of the injector. Since the time lag ct depends upon the chamber 
pressure it varies with time in the same manner as the chamber 


pressure does. In fact Tsien (4) has shown that 


Oar BG =e ea ee) 
——— ct ry tS 
dt e 
miere _ dlog HP Fo) 
r) lag 6 
(3) 


Now, as shown by Crocco’ ’, the continuity of mass flow 


through the combustion chamber may be written as 


3 
= + (Vetere) + eg (ze —3) —IMe > oy= Oo (1) 
where > is the reduced time lag a, and Z= a ° 
9, Og 


Furthermore, defining’ 
P= aoe ae P cn 
ep ZA pA Gs, 
the dynamic equilibrium of the propellant supply system is given by 


a 
dz 


=) + 





gr gti| y+ Pb sc (2) 


The presence of $) and w_ in Equations (1) and (2) indicates the 


coupling between the chamber and feed systems. 





To analyze the stability of the system the following solutions 


may be tried: 


Sz (S—s)2 Sz 
Pisye  Pcs-zB=Gesje , Mitei= mt)e 


D(z) 


where S=A+iw . Using these trial solutions and eliminating 


Se 


the common factor e Equations (1) and (2) represent a solution 


of the system if the two homogeneous equations below are simul- 


taneously satisfied: 


~33 ~3d 
(Sti-n+ne 396s) - e MCS) 


(3) 


I 
6) 


— 


Pcs) + | Js +3 0P+t)+0] m (5) 


it 
Q 


(4) 


For non-trivial solutions the determinant of the coefficients of 
Equations (3) and (4) must vanish, that is 


-s$ -sé 


S+i-n+nNe eee = 
2G60)so0 ©) 
te Js+ = GP =i +} 

This is a transcendental equation in the complex variable Ss _, the 
roots of which determine the stability of pressure oscillations. If 
any root of G’(s)=o possesses a positive real part the system 
is unstable. The Nyquist criterion may be applied to GG) to 
determine the existence of such roots (5) In this method the variable 
Ss traces a contour enclosing the right-half complex plane. As 5s 
moves in the manner shown in Figure 2, thatis, along the imaginary 
axis and a large half circle to the right, the behavior of <G/(s) is 


plotted and the number of complete revolutions noted. This number 


-7Je- 


is the difference between the number of roots and poles of Cs) in 
the right hali s-plane. If the number of poles can be determined 
independently the number of roots is known and the stability determined. 
The term ie in G'<s) makes use of the Nyquist diagram 

difficult. A complete investigation would reauire plotting the function 
fm <3) for many values of § , the reduced time lag. A different 
technique, proposed by M. Satche (6) and developed by Tsien ‘*) ; 
requires the function G‘(s) to be separated into two parts, the first 
consisting only of ts and the second of the remaining terms. In 


this case G’(s) may be expressed as a combination of two deter- 


minants 
S+1-/n ro) 
Os) = | (6) 
P JIS + ae at 
ry zal | 
Ds) | | (7) 
P wee Brees =) 


so that Equation (5) may be written 


—S& 
‘ho ee Di (s) = ° (8) 
-s% hie. 


A 
F ollowing Tsien' *) incalling e = eGo) and 





Equation (8) becomes 
See ie AC aaa) (9) 
and 9,(5) and G2(s) may be plotted separately. For neutral stability 


R2(s)=0O and Gy ao) becomes a unit circle on the Satche diagram; 


all unstable roots lie within the unit circle. Assuming that G(s) 


has no poles in the right-half plane then if 96.:s) either passes 
through the unit circle or encircles the origin the system is 
unstable, while the converse is true. The function G(s), however, 
may have poles in the right-half plane because of possible zeros in 
the denominator determinant 9,(5) . This number of poles may be 
found by applying the Nyquist criterion to 1 (5s) and this number 
must be added to whatever number of turns 9,(5) makes about the 
origin. 

As an example of the above criteria applied to a specific case 
a motor with the following parameters may be considered: 

A= 1.0 J «= 148 m2 0.6 P=1.0 


Expressing 9,(s) and J9.,(s) as follows: 


Gi ae + Ae + ae Sy (10) 
Do os) : aes a are (11) 
the coefficients become: 
dis [i+ 5(P +g) ]0-m d= Panfis 4 (P+ 4)] 
d= [ Jci-n) +14 = (P= i) Que a) (12) 
A CaN J 


For this particular motor 


Co) Co) 
on i = 25 
(1) 
= aq 
d, 3.) d, = 90,9 
(2) 
Gd. =15 


The characteristic equation for s is 


_ 55 (eues 24 |. Sia 
e. ae — = Oo 


e.5°26 O19 Ss 


> 





For the case where S='w , along the imaginary axis, 


: ee 
L+ 3B.) eu — 1.5 ay 
Q, Crew) ae 
2.5 =) Gronuee 


O.946w ~ z.5 i335 me 6a5 . 
= — ee Voc 


S254 6.8!m@ 62S >> o eras 


The Satche diagram for this function is shown on Figure 3. The 
function 9,(iu) crosses the unit circle and indicates instability 
for values of w below approximately 0.9. For large values of 
s off the imaginary axis, let S=R ae - Then o,.(s) behaves 
as - eS andas 9 decreases from 1/2 to -w/2 the 
curve progresses clockwise in a large arc from -o% to +~1 
in the left-half plane without encircling the origin indicating no 
further instability. 

References 4 and 5 contain petilincalia a stable and unstable 


motor behavior as determined by analyses similar to the above. 


- 9 ¢ 


il, FEEDBACK STABILIZATION 


The unstable operation of a motor with fixed parameters may be 
corrected by means of a simple feedback servomechanism, consisting 
of a combustion chamber pressure pickup instrument, sampling circuit, 
amplifier, servomotor and a variable capacitance. A feedback system 
in a stability analysis provides a mathematical relation between the 
combustion chamber pressure and the variable capacitance and intro- 
duces a new coupling term in the relations for the simple circuit with- 
out feedback. The relation between the capacitance « (z) and the 


pressure may be written symbolically as 


d 3 
in = —_— a 
eS ia gee i LF $03) 


where F_ S and F, = are linear differential-integral operators 
of the amplifier and sampling circuit, respectively. The above rela- 
tion indicates that the line capacitance changes with time in accordance 
with the amplifier signals whose action is controlled by a sampling 
circuit, which interprets the combustion chamber pressure impulses 


in a certain manner. 


A new set of simultaneous differential equations arises from the 





introduction of the feedback loop §9), 
d 
- + (1-n) > ony ( 2 > Sy - Ww (2-4) = 
z 
ae Busan ska , a 1 AKCRD 
 ° | (P+ = | fe + Po + JERS mae. ret 
c 9 \F < (2 | K (2) = 0 
d ' dz 


See 


Using trial solutions 


St sé st 
Kee) =n (s)veme Space) CP CS) Semmens = eo) 


a set of homogeneous simultaneous equations similar to (3) and (4) is 
obtained, 
sd -—s& 
(Sti-n+ne )$(s)-e MCS) =o 
\ \ 2 \ 

Pcs) + Js sf ee | mics) + lus + 4 (P+ 4)s| KX(s)=0 

F.(s1 F, (5)$ cs) —- K (5) = 0 
where F,(s) and F,(S) are the transfer functions of the sampling 
circuit and the amplifier circuit, respectively. The stability of the 


system now depends on the roots of the third order determinant, 


formed from the coefficients of the above set of equations 


-55 —~s$ 
S+l-n+ne -e © 
P Wee (era US (ae (13) 
F,(S) Fz¢s) ° = 


This determinant may be conveniently expressed in terms of 1G) 


Equation (6), D.(5) Equation (7), and a third minor defined as 2D¢ (s) 


| Oo 
One (14) 
Js+ os eee | Js + (P44) 


The new characteristic equation obtained by setting the stability 


determinant (13) equal to zero is 


~sd 
min e : [Dr +e, sits) Deis) | = Oo (15). 


For stability analysis using the Satche diagram Equation (15) may 
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be written 


- $2 D, cs) 

ae J,6s) + Fics) Fe(s) Dycsy 1 * 
where 9,(5) and 9,(5) now correspond to the first and second 
terms, respectively, of the above equation and may be plotted 
separately on the Satche diagram. Stable motor operation is to be 
assured by choosing the form of *,cs) and F,(s) so that the 
stability criteria are satisfied. It should be noted that D,¢s) , Dis) 
and {.:s) are functions of the fixed motor parameters and are 
independent of the feedback functions. 

Generally, instability in the basic motor arises from the inter- 
section of the 9,('w) curve with the unit circle as in Figure 3, 
where the reduced angular frequency is in the order of unity at the 
point of intersection. If for higher frequencies the function 9. 6) 
shows no further instability as is generally the case it is necessary 
only to modify 9,(iw) to move clear of the unit circle for small 


values of « and to cause the gain of the amplifier to disappear at 


higher frequencies. 








IV. SAMPLING CIRCUIT FOR FEEDBACK SYSTEM 


The functions of the feedback circuit in this investigation will 
be divided into two parts: one, to provide stabilization and control 
the gain at different frequencies and the other, to transform the 
pressure signal received from the combustion chamber in such 2 
manner that the compression stroke of the variable-capacitance will 
be of the same frequency and the return stroke will take considerably 
longer time. A sampling circuit will be responsible for carrying out 
this latter task while ‘is stabilizing function will be carried out by the 
amplifier. 

As a representative example, let the pressure signals received 
by the sampling circuit from the combustion chamber be sinusoidal, 
and the signals which the sarnpling circuit delivers to the amplifier 
be represented by an asymmetric function, as follows: 


Into the sampling circuit 


sin z 


Out of the sampling circuit 


(cla sed interval] QO< 2<7 
f (2) 


é-T 
<q- i 


= — cos 





interval W< 2C2@ng 


The input and output functions are shown graphically on Figure 4. 
Over a total interval of time 2% q there are q pressure cycles in 
the combustion chamber for each cycle of motion in the variable 


capacitance. 


- eee 


The function 42) may be expanded in a Fourier series 


Le 8) 
i a's iin 46 
Hess - y Ga cos q t+ b sin am 
k= 
The coefficients are 
arg 
al ee | fez) cos K z dz 
WUQ Jo 2 
21rq : 
a\ 5 — | $(z) Eo ¢ dz 


TS a 


Replacing f{z) by the required output functions in the above equations the 


following general expressions for the coefficients are obtained 














| = 
pee.” commie: Joe | scas ==" eg. oe 
MUS) " S| TG 2o-9 4 
o 
: | 
= K Sim k Th ~ ee ae 
qT a q+ ke Qt _ 
CZa as 
T cig 
\ ; 
b= a COS @ sim 2 ae = 66s =" sin eee 
a ee 2 ze) 2q-1 Q 
k kk 
i+ cOos—T aa 
: A q COs q 
“<r q* a L.* 2 oa 
(zQ- Ge a ~ 


The second expressions for a, and 2), above may be used to 


evaluate the coefficients when q7#k For the case where ee 


the integra! expressions used directly are more convenient. 





Anticipating the use of an illustrative exarmple the coefficients 
for the first eight terms of the Fourier series using a value of q=4 


have been calculated and are tabulated below. 











0. 350 





0.821 





0.226 “0.1731 






0.174 0.0643 






0.125 






0 





0.0794 ~0.144 













0. 0420 -0.0535 





0.0154 “0.0303 





0 0.027 





Now, it is necessary to obtain the transfer function of this 
sampling circuit in terms of the complex variable s eventually to 
investigate the stability of the syatern. The transfer function is the 
ratio of the Laplace transforms of the output and input functions of 
the sampling circuit. 


The Laplace transform of the output function is 


oe] Co 


Sz ‘ —Sz 
4(Z)e dz = > (a,cosk2+ on Sram e dz 
° o k=t “ 
co oo MK °° Q 

5 q 
as 5-Vorers! " es ~ ) Ses gs ee + Bk) 
eee ee Ge 
Similarly for the input circuit 
os 
; —Sé 
Sine © Slee —— 
5S + 4 


@ ) 
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Therefore, the transfer function of the sampling circuit may 


be written as 


oo aN ; 
= ee ES oo 
FCs) ) “EA we, (a,qs cK) (17) 
Kz q* 


Introducing aterm CcsS in (17) to take into account the damp- 
ing which exists in a real circuit because of resistance elements the 


function F, (S) will be modified as follows: 
Rf ees S° 


“Crane. 
k=, 4 is ie Sess coe 


Fis) modified = (A.qs + bs) (18) 


It is now posstble to proceed to the development of the appro- 
priate amplifier transfer function and to determine the overall 


stability of the system. 
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V. AMPLIFIER CHARACTERISTICS AND STABILITY ANALYSIS 


An appropriate procedure for determining the amplifier 
characteristics is (5), (1) to find the approximate behavior of F,(s) 
for S=iw<cc. suchthat 9,((w) does not intersect the unit circle 
and (2) to determine the function F,:5) so that it has the prescribed 
behavior for small values of ~ and vanishfor large w . The 
first part of the procedure may be carried out in the following 


manner. Each of the terms in the expression for 9.¢°), 


D , Cs) 
9,093 = = —— (19) 
OD, (s) + FCS) F2¢s) O¢ Gs) 


obtained frorn Equation (16) may be expressed as a power series in 


Thus 
SiGe d "+ dos + de 
Dos) a aes o cee 
and eo 2.5 + Boe (20) 
4 \ 
where ae (ess 2) 2 Bene (21) 


The B, are obtained by expanding the determinant of Equation (14). 


Theterms Js) D,(s) and Dy (s) depend only on the fixed 


parameters of the basic motor. Equation (18) representing the 
transfer function fF, (s) of the sampling circuit may be expanded 


as follows 





On| 
ea qb, (+ 20545) (1+ <— s) 
oe ee 








24, 
K=t Kk }+ = es ae s*) 
wn 194 
vad qb, 2 a . 2,2 eo 
= > (i+2cS5+S°)|(——(2cs+s°)+ — (2ecs4 8°) 4... (. 5) 
kk * Vea b, k 


Kaz 














3 aq\| q ie 
a 2 fo) 3 4c —| Se (22) 





It is appropriate to choose the following form for the amplifier 


transfer function 


me (1), Co) (1) 
Es) = {_ es =F a Se se S (23) 


where the {/' are unknown constants. The function 9.5) may 


now be written 
Gs ache ees 


{-«) 


cor Gs) 2 £ (oy 1) 2 (24) 
d, +d, s + (A,+A,5+A,5 ) ( > + f, +4. s)(B,s+ B.S ) 


9,6) z 


aS 
where, for brevity, A, represent the coefficients of s in the approxi- 
mate expression for F,(s)  , Equation (22). The reason for choosing 
.the form for F2(s) of Equation (23) becomes clear from inspection 
of Equation (24), since, for s = 0 


S a 


UO) 7 eee emer rere ie ee 
: Ane ies = 


and 9g,(0) remains finite. Now the requirement that the 
curve lie outside the unit circle for small WG may be met as 


follows. Prescribe 9,(iw) by the polynomial 
9, (iw) = So + Cuca Vo Caucuses (25) 
choosing the %, so that stability will be obtained. Then there are 


aw 








. 19 - 


two expressions for 9,(5) validfor small s | namely 


2 


G.05)= % + %,3 + 6,5 (26) 
and Equation (24). Since these two expressions are power series 
expansions in s they are equal term by term. The results of setting 
coefficients of like powers of s equal to each other are cipvinviend 
relating the unknown constants iy and the known coefficients from 
which the values of ie may be calculated. Thus, after performing 


the necessary operations and retaining terms in s of no higher 


power than a the following expressions for the {- are obtained: 











(-1) 1 Sh, Co) 
+. ‘eos (-<: aa d, (27) 
ee K 
(0) ou 
(0) | el v, qb, (0) 
i = ~ 2 av, -2- = (a) eo: (28) 
nee . ; zt 
k=1 K 
db, , ©) q Kd qb, . G) (1) 
oJ fe(-SpS8)-0b BO 8 
key - Ke 
(rv) CO a a 2 1 
a = of ox a qb. ss ea ~*~ (2c pe: ) = 5 [1-4 2 | 
yoke Li —_ : O. as ele : 
K : 
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The second part of the procedure for determining the complete 
amplifier transfer function is to obtain an expression which agrees 
with Equation (23) for small values of s and will vanish for large values 
of 3 so as to leave the original Satche diagram unchanged for large s. 
Consider the following choice of [f,(5) 


(-9) Zz 
fie as+b 

Fts)= 2 (+ =~\ (30) 
S 


{+ hs 








where a, b, and h are constants. For small values of s 


£ (-1) 
me Z 
S 


eras + 9b-: ) 





ese 


and the coefficients a and b must become 


(o) (4) 
{os te 





ic ro ae co pa) 
in order that -,(s) may agree with Equation (23). For large s, 
i 
eecS) will behave as ae , causing the second term in the 
S 


denominator of 9,(s) in Equation (19), that is, the product 
feo) F, (5) D6 (s) ,» to behave as s° so that the value of Q,¢5) 
becomes nearly that of the basic motor without feedback. The value 
of h must be such that stability will be obtained for all values of s 
and will be determined largely by trial and error. 

The exact expression for the function 9,(5) becomes then, 


for all values of s 


) 
— ( ae fer, ae ese) 
G,(S)= = = 
(i+2e545}(go,s+kb,) £5" (14 054 bs") 
ep) @) 2 = 
d, +0, $ a Zz 4 (B,s+ Bis) 
Ko + 2cS4 Ss" S>  oi+hs 


oie (32) 
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To illustrate the application of the stabilization procedure 
consider the unstable motor for which the Satche diagram is shown 
on Figure 3. To provide stabilization at low values of ~ the curve 
of 9,(iw) must be moved away from unit circle. Choosing the 


following values of x. 


\ 


8, = °2. 25 | v,2 22.0 <¢,2 -0.25 


and using Equation (25) for 09, (iuw) the dashed line on Figure 3 is 
obtained. From Equations (27), (28) and (29), using a value of 

c = 0,1, and summing over the first eight terms of the Fourier 
coefficients, the { . are 


(- 


T, 


(0) 


"20.304 9 f° 0.032 fs 5.34 | 


which fixes the feedback loop for small values of w . To obtain the 
general expression for the amplifier transfer function the constants a 


and b are, from Equations (3}), 
a = 3.395 b= 17.48 


and F., cs) becomes 


OZ ow Cees oes Weersice | 


a 


(ey\ Ss 
- > L+H Ss 


Choosing arbitrarily a value of h=1 the function 9,5) becomes 


Se ay Ae 


OQ, (S)= 


25a Vs wee eS. —— 
4 2 
—_ Te Ss =, 
= a ti SsS-+ 


& 
(\+0.25+S5 \(4aps+ kb, ) 0304 [t+33955 + 17.488" 
i+ Se 


Jus + 1.55") 


Yor the condition where S=tw 


=—(\— Ses 45 seme co } 


a0 w )= S 
(i~ iow OZiw\(kb+ 4a,iw) \.034 0.304 ~5.32w~ 
Z.540diw+ — > 


ae SS aay pe 


c 2 
= —- WwW + Q,.ckw 


meor Us(i+ wy 
kza iG 


This function is plotted on Figure 5 and shows failure to provide 
complete stability since it cuts through the unit circle. 

For a value of hW=10 however, the function 9,(:w) stays 
clear of the unit circle for small values of q@w #, as shown on 
Figure 6 (on this plot also are shown the unstable curves of the motor 
without feedback and the curve for the approximate behavior of 
with feedback for small values of w ). As w increases, 9. iw) 
moves down toward —co. ; as s traces the arc in the right-half 
plane shown on Figure 2, 9,(5) proceeds in a large clockwise arc 
in the left half-plane and returns to the positive imaginary axis and 
finally traces a mirror image of the contour below the real axis. The 
closed contour, therefore, does not enclose the unit circle. 

To check for unconditional stability the Nyquist criterion must 
be applied to the denominator of 9,(5) , namely D, 0s) + Fc) Ges) Deosy, 
4.8 shown on Figure 7, the plot of this function makes no revolutions 
around the origin demonstrating unconditional stability for all values 
of time lag. 

Of special interest is the behavior of 9, (iw) in the vicinity of 
o> =0.9. For frequencies close to this yalue the proximity of the 
curve to the unit circle shows that pressure oscillations in the com- 


bustion chamber may be poorly damped. 








= a5 °- 
VI. CONCLUDING REMARKS 


The scheme of introducing a sampling circuit in the feedback 
loop to transform the pressure signals in a certain manner and 
alleviate probable mechanical difficulties in the variable capacitance 
does not appear to make the problem of designing a stabilizing 
circuit more difficult than usual. The question now arises as to the 
possibility of designing a practical sampling circuit with real 
components that will transform symmetric signals to non-symmetric 
signals. It is necessary to introduce an actual nonlinear circuit into 
the problem and to analyze its stability and performance. In all 


probability the technique of analysis must be modified. 
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